The globular cluster NGC 1978 in the Large Magellanic Cloud 



Alessio Mucciarelli 

Dipartimento di Astronomia, Universitd degli Studi di Bologna, Via Ranzani, 1 - 40127 

Bologna, ITALY 

alessio .mucciarelli@studio .unibo . it 

Francesco R. Ferraro 

Dipartimento di Astronomia, Universitd degli Studi di Bologna, Via Ranzani, 1 - 40127 

Bologna, ITALY 

f rancesco . f erraro3@unibo . it 

Livia Origlia 

INAF - Osservatorio Astronomico di Bologna, Via Ranzani, 1 - 40127 Bologna, ITALY 

livia. origliaObo . astro . it 
Flavio Fusi Pecci 

INAF - Osservatorio Astronomico di Bologna, Via Ranzani, 1 - 40127 Bologna, ITALY 

flavio . f usipecci@oabo . inaf . it 

ABSTRACT 

We have used deep high-resolution Hubble Space Telescope ACS observations 
to image the cluster NGC 1978 in the Large Magellanic Cloud. This high-quality 
photometric data set allowed us to confirm the high ellipticity (e ~0.30±0.02) 
of this stellar system. The derived Color Magnitude Diagram allowed a detailed 
study of the main evolutionary sequences: in particular we have detected the so 
called Bump along the Red Giant Branch (at V555=19.10±0.10). This is the first 
detection of this feature in an intermediate-age cluster. 

Moreover the morphology of the evolutionary sequence and their population ra- 
tios have been compared with the expectations of different theoretical models 
(namely BaSTI, PEL and Padua) in order to quantify the effect of connec- 
tive overshooting. The best agreement (both in terms of morphology and star 
counts) has been found the PEL (Pisa Evolutionary Library) isochrone with 
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Z=0.008 (consistenly with the most recent determination of the cluster metallic- 
ity, [M/H]=-0.37 dex) and a mild overshooting efficiency (A os =0.1). By adopting 
this theoretical set an age of r=1.9±0.1 Gyr has been obtained. 

Subject headings: Magellanic Clouds — globular clusters: individual (NGC 1978) 
— techniques: photometry 



1. Introduction 



The Large Magellanic Cloud (LMC) is the nearest galaxy of the Local Group with 
evident star formation activity. The LMC has a very popolous system of globular c l usters 
that covers a large range of metallicities and ages (see lOlszewski. Suntzeff fc Mated (119961 ) 
and references therein). We can distinguish an old population, coeval with the Galactic 
globular cluster (GGC) system, an intermediate (1-3 Gyr) and a young population, ranging 
from a few million years to 1 Gyr. Hence the stellar populations of the LMC cluster system 
represent an important laboratory to study and test stellar evolution models in a different 
age-regime with respect to the Galactic system. 

The first s urveys devoted to stu dy the properties of these clusters have been ca rried out 



in the optical (Ivan den Berghl Il98ll ) and in the near- infrared (IPersson et al.l 119831) spectral 



range and were based on measurements of their integrated colours (jSearle. Wilkinson, fc Bagnuolo 



19801 : lElson fc Falll 1 1 98 5f ) . However, only in the last decade the main properties of the re- 



tion photometry of indivi d ual stars JVallenari et al. 
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20061 ). These works have produced a number of interes t ing results: (i) the oldest clus 
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the LMC are coeval wit h GGCs (ITesta et al.l 11995c iBrocato et al.lll996l : lOlsen et al. 



The recent census of the old LMC globulars presented by 
a total of 15 objects with ages > 10 Gyr; (if) the dominant 
population is constituted of ~100 young and intermediate age clusters, with ages <3 Gyr; 
(Hi) the absence of objects in the huge age range from ~3 to ~ 13 Gyr (the so-called Age 
Gap): The Age Gap problem has been l argely investigated 



clust ers in the Age Gap actively searched (IGeisler et al. 



1997 



'see i.e. 


Be 


cki et 


al. 


2004 


) and 


Rich et al. 


2001 


Mackev et al. 



20061 ). however up to date, only one candidate, namely ESO 121-SC03, possibly falls in such 
an age range. In spite of these efforts, the overall formation history of the LMC is still 
an open issue, since a firm age-metallicity relation (AMR) has not been established, yet. 
Indeed, homogeneous and accurate metallicities and ages for a significant sample of LMC 
clusters are mandatory to derive a reliable AMR. 
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Accurate ages from the measurement of the Main-Sequence (MS) Turn-Off (TO) region 
are still sparse and very model (i.e. isochrones) dep endent. The only homoge neous age-scale 



available still relies on the so called s-parameter (jElson &: Falll Il985l . Il988l ). an empirical 



quantity related to the position of the cluster in the dereddened (U-B) vs (B- V) color-color 



diagram. This paramete 
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Most chemical abundance determinations for LMC clusters are based on low-resolution 



spectra ({Olszewski et al. 



19911). integrated i nfrared spectroscop y (jOliva fc Origliaj Il998[ ) or 



photometric techniques (IDirsch et al.l |2000| ; lLarsen et al.l 120001) . Detailed chemical abun- 



dances from high-r esolution spectra are still sparse (IHill et al.l |2000| ; iJohnson et al 
Ferraro et al.ll2006l ). 



2006 



With the ultimate goal of constructing a homogeneous age-metallicity scale for the 
LMC clusters, we started a program which makes use of the last generation of instruments 
(imager and multi-object spectrograph) in order to perform an appropriate study of stellar 
population, age, metal content and structural parameters for a number of pillar clusters. In 
this paper we present the results for one of the most massive and luminous globulars in the 
LMC, namely NGC 1978. This clust er has been conside red for years a peculiar object (i) 
because of its high-ellipticity (e=0.3, Fischer et al. 19921 ). significantly larger than the one 



typically measured in stellar clusters ( Goodwinll997 ). and (ii ) because it has been suspected 



to ha rbour a che mically inhom o geneo us stellar population (lAlcaino et al.l Il999l ; IHill et al. 
20001 ) . However, iFerraro et al.l (120061 ) has recently presented high resolution spectra for 



eleven giants in this cluster, obtaining [Fe/H]=-0.38 dex, with a very low dispersion (<r=0.07 
dex), firmly excluding the presence of a significant metallicity spread in this cluster. Here 
we present a high precision CMD of the cluster based on observations obtained with the 
Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope (HST). The 
morphology and the population ratios along the main evolutionary sequences have been 
used to quantify the effect of convective overshooting. Finally, by co mbining the precise 
measure of the MS-TO with accurate estimate of the cluster metallicity (IFerraro et al.ll2006l ) 
we also provide an accurate estimate of the cluster age. 



2. Observations and data analysis 



2.1. Observations 



The photometric dataset consists of high-resolution images obtained with ACS@HST 
(300 sec and 200 sec long) through the F555W and F814W filters. These images have 
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been retrieved from ESO/ST-ECF Science Archive (Proposal ID 9891, Cycle 12). The 
observations have been obtained with the Wide Field Channel (WFC) that provides a field 
of view of ^200" x200" with a plate scale of 0.05 arcsec/pixel. The WFC is a mosaic of 
two CCDs, both with 4096x2048 pixels separated by a gap of ~50 pixels. The first chip 
has been centered on the cluster center, while the second chip sampled a contiguos field. 
All images were reduced with the ACS/WFC pipeline, in order to perform bias and dark 
subtractions and fl atfield correct ion. The photometric reduction was performed by using 
the DA OP HO T-II JstetsorJ [l987|) Point Spread Function (PSF) fitting method. The output 
catalog includes about 40,000 objects from the first chip and more than 10,000 from the 
second chip and it has been calibr ated in the ACS/WFC Vega-mag system, following the 
prescriptions of iBedin et al.l (120051 ) . Finally, the ACS catalog has been astrometrized in 



the 2M ASS astrometric system by matching the IR catalog presented in iMucciarelli et al. 
f|2006h . 



2.2. The Color-Magnitude Diagram 

Fig. 1 shows the calibrated color-magnitude diagram (CMD) for the chip centered on 
the cluster. Stars in the brightest portion of the Giant Branches could be saturated and/or 
in the regime of non linearity of the CMD. Hence for stars brighter than Vs55=17.6 (this 
magnitude level is marked with an horizontal dashed line in Fig. 1), magnitudes, colours 
and level of incompleteness are not safely measured. This CMD (reaching the magnitude 
limit of V555 ~26) shows the typical evolutionary features of an intermediate-age stellar 
population, namely: 

(1) the brightest portion of the MS at V555 <21 shows a hook-like shape, typical of the 
evolution of intermediate-mass stars ( M > 1.2M J which develop a con nective core 0. In 



particular, the so-called overall contraction phase (ISalaris fc Cassisil 120061 ) is clearly visible 
between the brightest portion of the MS and the beginning of the Sub-Giant Branch (SGB) 
at V 555 ~20.9. 

(2) the SGB is a narrow, well-defined sequence at V555 ~20.7, with a large extension in 
color (<5(V555-/si4) ~0.6 mag). The blue edge of the SGB is broad and probably affected by 
blending, especially in the most internal region of the cluster. 



1 Note that the width of the distribution in color of the bright portion of the MS (a^y-i) ^0.05 mag) 
turns out to be fully consistent with the observational errors estimated from the completeness experiments 
(cry ~ a j ~0.03 mag, corresponding to <J(y-i) ^0.04 mag). 
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(3) the Red Giant Branch (RGB) is fully populated; this is not surpris ing since this 
cluste r has already experienced the RGB Phase Transition (see the discussion in lFerraro et al. 
r |2004h and lMucciarelli et all (bood )). 



(4) the He-Clump is located at V555 ~19.1 and (V555 — I&u) ~1.15. 

Fig. 2 shows the CMD of the external part of the ACS@HST field of view (corresponding 
to r>140" from the cluster center). This CMD can be assumed as representative of the field 
population surrounding the cluster. In particular, the CMD shows two main components: 

(1) a blue sequence extended up to V555 ~17. 

(2) a SGB which merges into the MS at V555 ~22.2, corresponding to a population of ~5 
Gyr. We interpret this feature as a signature of the major star-formation episode occurred 
5-6 Gyr ago, when LMC and Small Magellanic Cloud (SMC) were gravitationally bounded 
(feekki et al1l2004h . 



2.3. Completeness 

In order to quantify the degree of co mpleteness o f the final photometric catalog, we used 



the well-know artificial star technique (jMated Il988l ). and simulated a population of stars 



in the same magnitude range covered by the observed CMD (excluding stars brighter than 
^555=17.6, corresponding to the saturation level) and with (V555 — Jgu) ~0.8 mean color. The 
artificial stars have been added to the original images and the entire data reduction procedure 
has been repeated using the enriched images. The number of artificial stars simulated in each 
run (~ 2,000) are a small percentage (~5%) of the detected stars, hence they cannot alter the 
original crowding conditions. A total of ~250 runs were performed and more than 500,000 
stars have been simulated. In order to minimize the effect of incompleteness correction, we 
have excluded the very inner region of the cluster (r<20", where the crowding conditions are 
most severe) from our analysis. In Fig. 3 the completeness factor = (defined as the 
fraction of recovered stars over the total simulated ones) is plotted as a function of the V555 
magnitude in two different radial regions, namely between 20" and 60" and at r>60" from 
the cluster center, respectively. 



2.4. The RGB-Bump 



The extended and populated RGB in NGC 1978 gives the possibility to search for the so- 
called RGB-Bump. This is the major evolutionary feature along the RGB. It flags the point 
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when the H-burning shell reaches the discontinuity in the H-abundance profile left by the in- 
ner pene tration of t he convection. This feature has been predicted since the early theoretical 
model s (llbenl 119681) bu t observed for the first time in a globular cluster almost two decades 
later (IKing et al .11 1 9 84f) . Sinc e that first detection the RGB-Bump w as identified in several 



GGCs (IFusi Pecci et al.lll990l : iFerraro et all 



1999 



in the Loca l Group (Sextant (jBellazzini et al.l 120011 ). Ursa Minor ( 



Zoccali et al.lll999) and in a few galaxie s 



Bellazzini et al. 2002). 



Sagitt arius ( Monaco et al. 20021 )). Accordingly with the prescriptions of Fusi Pecci et al. 
(119901 ). we have used the differential and integrated luminosity function (LF) to identify the 
magnitude level of the RGB-Bump in NGC 1978. In doing this, we have (1) selected stars 
belonging to the brightest (V555 <20.6) portion of the RGB; (2) carefully excluded the bulk 
of the He-Clump and AGB stars by eye; (3) defined the fiducial ridge line for the RGB, 
rejecting those stars lying at more than 2a from the ridge line. Fig. 4 shows the final RGB 
sample (more than 600 stars) and both the differential and integrated LFs. The RGB-Bump 
appears in the differential LF as a well defined peak at V 5 ^ p = 19.10 and it is confirmed in 
the integrated LF as a evident change in the slope. 

For both LFs the assumed bin-size is 0.1 mag; in order to check the uncertainty in the 
Bump magnitude level, we have tested the position of this feature by using LFs computed 
with different binning. The impact of the selected bin-size is not crucial: a difference of 0.2 
mag corresponds to a variation <0.05 mag in the detection of RGB Bump. By considering 
the intrinsecal width of the peak in differential LF, we estimate a conservative error <0.10 
mag. 

Finally, we note that the RGB Bump is brighter and reddest than of the bulk of the 
He-Clump and the latter merges into the RGB at faintest magnitude (V555 ~19.3, see Fig. 
4; hence the possibility of contamination is negligible. 



3. The cluster ellipticity 



Most globu lar clusters in the Gal axy show a nearly spherical shape, with a mean ellip- 
ticity I e=0.07 Jwhite fc Shawllll987h and more than 60% with e <0.10. One of the most 
remarkable exception is represented by u Centauri that is clearly more elliptical than the 
other GGCs: its ellipticity is 6=0.15 in the exte rnal regions with a evident decrease in the 
inner regions, with e=0.08 (jPancino et al.l 120031 ) . Conversely, the LMC cl usters (as well as 
those in the SMC) show a stronger departure from the spherical simmetry. iGeisler fc Hodge 



2 Note that ellipticity is defined here as e=l-(b/a), where a and b represent major and minor axis of the 
ellipse, respectively. 
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( 11980l ) estimated the e llipticities of 25 popolous LMC clusters, finding a mean value of 
e=0.22; iGoodwinl Ji997m obtained a lower average value for the LMC clusters (e=0.14), but 
still higher than the mean ellipticity of the GGCs. Moreover, in the LMC the presence of 
many double or triple globular clusters has been interpretated as a clue of the possibility 
of merger e pisodes between su bclusters with the result to create stellar clusters with high 
ellipticities ( jBhatia et al.lll99lf). 



Previous determinations (IGeisler fc Hodgdll980l ; iFischer et al.lll992l ) suggested large values 
of ellipticity for NGC 1978. 

We have used the ACS catalog to derive a new measurement of the ellipticity of the 
cluster, in doing this we computed iso de nsity curve s and a dopting an adaptive kernel tech- 
nique, accordingly to the prescription of iFukunagal (ll972l ). In doing this we have adopted 
the center of gravity of the cluster computed using the near-infrared photometry obtained 
with SOFI (jMucciareUi et al.ll2006l ). The isodensity curves have been computed using all the 
stars in the first chip with V555 <22 (approximately two magnitudes below the TO region) in 
order to minimize the incompleteness effects §. Finally, we have fitted the isodensity curves 
with ellipses. Fig. 5 shows the cluster map with the isodensity contours (upper panel), 
the corresponding best fit ellipses (central panel) and their ellipticity as a function of the 
semi-major axis in arcsecond (lower panel). No evidence of subclustering or double nucleus 
is found. The average value of the ellipticity results e=0.30 (with a root mean square of 
0.02), without any radial trend. This value is in good agreement with the previous estimates 
(IGeisler fc Hodgdll980l ; IFischer et al.lll992l ) and confirms the surprisingly high ellipticity of 
NGC 1978. 



4. The cluster age 



The determination of the age of a stellar population requires an accurate measure of 
the MS TO and the knowledge of the distance modulus, reddening and overall metallic- 
i ty. For NGC 1978 we used the recent accurate determination of [Fe/H]=-0.38±0.02 dex 
(IFerraro et al.l 120061 ) and < [a/Fe] > almost solar (Mucciarelli et al., in preparation), based 
on high-resolution spectr a, to derive th e overa ll metallicity [M/H]. In doing this, we adopted 
the relation presented by ISalaris et al.l (119931 ): 



[M/H] ~ [Fe/H] + log (0.694 • 10 <[a/Fe]> + 0.306), 



3 Note that a different assumption on the magnitude threshold does not affect the result. 
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obtaining [M/H]~ -0.37 dex. 

In the case of intermediate age stellar systems, the measurements of the age is complicated 
by the presence of a convective core, whose size needs to be parametrized (A os ) 

B We 

then use different sets of theoretical isochrones with different input physics, in order to 
study the impact of the convective overshooting in reproducing the morphology of the main 
evolutionary sequences in the CMD. 



BaSTI mo dels: BaSTI ( A Bag of Stellar Tracks and Isochrones) evolutionary code 
described in iPietrinferni et al.l (120041 ) computes isochrones with and without the inclu- 
sion of overshooting. The overshoot efficiency depends on the stellar mass: (1) A os =0.2 
for masses larger than 1.7 M ; (2) A os = 0.25 ■ (M/M© - 0.9) for stars in the 1.1-1.7 
M range; (3) A os =0 for stars less massive than 1.1 M & . 



Pisa models : PEL (Pisa Evolutionary Library, ICastellani et al.l (120031 )) provides an 
homogeneous set of isochrones computed without overshooting and with two different 
values of A os , namely 0.1 and 0.25. 



Padua models : in these isochrones (IGirardi et al.ll2000l ) A^ dua =0 for stars less mas- 
sive than 1 Mq, where the core is fully radiative. The overshooting efficiency has been 
assumed to increase with stellar mass, according to the relation A^, ° = M/M Q — 1 
in the 1-1.5 M range; above 1.5 M a costant value of A^ dua =0.5 is assumed. Note 
that this value corresponds to A os ~ 0.25 in the other models, where the extension of 
the convective region (beyond the classical boundary of the Schwarzschild criterion) is 
measured with respect to the convective core border. 



4.1. The morphology of the evolutionary sequences 

From each set of theoretical models, we selected isochrones with Z=0.008 (corresponding 
to [M/H]=-0.40 dex), consistent with the overall metallicity of the cluster and we assumed a 



4 The overshooting efficiency is parametrized using the mixing length theory ( Bohm-Vitense 19581 ) with 



A os =l/ Hp (where H p is pressure scale height) that quantifies the overshoot distance above the Schwarzschild 
border in units of the pressure scale height. Some models as the "Padua ones" define this parameter as the 
overshoot distance across the Schwarzschild border, hence the A os values from different models are not always 
directly comparable. 
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distance modulus ( m — M)o ^18.5 (Ivan deu Berghlll998l ; IClementini et al.ll2003l ; lAlveal2004j ) 
and E(B-V)=0.10 (IPersson et al.lll983l ). However, in order to obtain the best fit to the ob- 
served sequences with each isochrone set, we left distance modulus and reddening to vary by 
< |10|% and < |30|% , respectively. Fig. 6 shows the best fit results for each isochrone set, 
while Table 1 lists the corresponding best fit values of age, reddening, distance modulus and 
the predicted magnitude level for the RGB-Bump. The best fit solution from each model 
set has been identified as the one matching the following features: (i) the He-Clump mag- 
nitude level, (ii) the magnitude difference between the He-Clump and the SGB and (Hi) the 
color extension of the SGB. The theoretical isochrones have been reported into the observa- 
tional plane by means of suitable transformations computed by usin g the code described in 



Origlia fc Leithererl (120001 ) and convolving the model atmospheres by lBessel. Castelli fc Plez 



(119981 ) with the ACS filter responses. In the following, we briefly discuss the comparison 
between the observed evolutionary features and theoretical predictions. 



• BaSTI models: By selecting canonical models from the BaSTI dataset, the best fit 
solution gives an age of 1.9 Gyr, with E(B-V)=0.09 and a distance modulus of 18.47. 
Despite of the good matching of the He-Clump and SGB magnitude level, and the RGB 
slope, this isochrone does not properly reproduce the shape of the TO region and the 
overall contraction phase (panel (a) of Fig.6). The best-fit solution from overshooting 
models gives an age of 3.2 Gyr, (m — M)q= 18.43 and E(B-V)=0.09, and matches 
the main loci of the evolutionary sequences in the CMD. In particular, this isochrone 
provides a better match to the hook-like region (between the MS and the SGB, see 
panel (b) of Fig.6). 



• PEL models: Panels (c), (d) and (e) of Fig. 6 show the best fit solutions obtained 
by selecting 3 different A os . In all cases (m — M)o=18.5 and E(B-V)=0.09 are used. As 
can be seen values of A os =0 and A OS =0.25 isochrones fail to fit the SGB extension and 
the hook-like feature, conversely a very good fit is obtained with a mild-overshooting 
(A os =0.1) and an age of r=1.9 Gyr. 



• Padua models: The best-fit solution gives r=2.2 Gyr, (m — M) =18.38 and E(B- 
V)=0.07 ( Panels (f) of Fig. 6). This isochrone well-reproduces the complex structure 
of the TO and the core contraction stage, as well as the SGB structure and the RGB 
slope. However, it requires distance modulus and reddening significantly lower than 
those generally adopted for the LMC. 
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From this comparison, it turns out that only models with overshooting are able to best fit 
the morphology of the main evolutionary sequences in the observed CMD. In particular, 
the best fit solutions have been obtained with the BaSTI overshooting model with r=3.2 
Gyr, the PEL mild-overshooting model (A os =0.1) and t—1.9 Gyr and the Padua model with 
\Padua= 25 (corresponding to A OS =0.25) and r=2.2 Gyr. 

However, it must be noted that none of these models satisfactorly can fit the observed 
Bump level, the BaSTI and Padua models being ?a0.1 and 0.3 mag fainter, respectively, and 
the PEL model £s0.2 mag brighter, perhaps suggesting that evolutionary tracks for stars with 
M>1M Q still need some fine tuning to properly reproduce the luminosity of this feature. 



4.2. Population ratios 

Since the comparison between the observed CMD and theoretical isochrones is somewhat 
qualitative, we also performed a quantitative comparison between theoretical and empirical 
population ratios: this yields a direct check of the evolutionary timescales. To do this, we 
define four boxes selecting the stellar population along the main evolutionary features in our 
CMD, namely the He-Clump, the SGB, the RGB (from the base up to V555 ~19.4) and finally 
the brightest (~ 1 mag) portion of the MS; these boxes are show in Fig. 7, overplotted to the 
cluster CMD. Star counts in each box have been corrected for incompleteness, by dividing 
the observed counts by the factor obtained from the procedure described in Sec. 2.3 (see 
also Fig. 3) for each bin of magnitude. 

Star counts have been also corrected for field contamination. To estimate the degree of 
contamination by foreground and background stars we have applied a statistical technique. 
We have used the CMD shown in Fig. 2 as representative of the field population. The 
number of stars counted in each box in the control field have been normalized to the cluster 
sampled area and finally subtracted from the cluster star counts. The final star counts per 
magnitude bin in each box have been estimated accordingly to the following formula: 

AT _ Nobs _ ^ 

-<> corr , 1 'i field j 

where N obs are the observed counts and Nf ietd the expected field star counts. We find 
iV M s=4331, N SG b=Q32, N RGB =450 and N Heci =3U, where N MS , N RGB , N SG b and N He _ cl 
are the number of stars in the box as sampling the MS, the RGB, the SGB and the He-Clump 
population, respectively. 

Uncertainties in the computed population ratios have been estimated using the following 
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formula 



0R 



VR 2 -vh + v 2 



N 



D 



where R = N/D is a given population ratio, N is the numerator and D the denominator 
of the ratio. The errors cr/v and <j£> for any population have been assumed to follow a 
Poisson statistics. In addition, in the error budget we also include the uncertainty due to the 
positioning of the box edges: note that a slightly different (± la) assumption in the definition 
of the box edge has little impact (typically 7-8%) on the star counts. This uncertainty has 
been quadratically added to the Poissonian error. 

On the basis of the boxes shown in Fig. 7 we defined four population ratios 0, as listed 
in Table 2: (i) N RGB /N SGB ; (ii) N RGB /N He _ cl ; (iii) N SGB /N He -ci] (iv) N MS / N {SGB+RGB) . 
For each selected model, corresponding theoretical population ratios have been estimated by 
convolving the isochrone set shown in Fig.4 with an Initial Mass Function (IMF), according 
with the prescriptions of IStraniero &: Chieffil (Il99ll ). In order to check the sensitivity of the 
popu lation ratios to the ad opted IMF, we have use d three diffe rent values for the IMF slope a: 
2.35 JSalpeterlll955T l. 2.30 fcroupalboOll ) and 3.5 JScalolll986h at M>1M . In the considered 
mass range (between 1 and 2 M ), the theoretical population ratios are poorly dependent 
on the assumed IMF, with a 16% maximum variation (between Scalo and Kroupa IMFs) for 
the Nms I N(sgb+rgb) ratios. Hence in the following the population ratios are computed by 
using a Salpeter IMF. 

As results, we found that BaSTI and PEL canonical models predict a lower (by <40%) of 
the N M s/N(sgb+rgb) and higher N RGB /N He _ C i and N SGB /N He _ C i (by <35% and <100%, 
respectively) population ratios with respect to the observed ones. Isochrones with high over- 
shooting (A OS =0.2-0.25) show an opposite trend, with higher (by <50%) Nms I ^(sgb+rgb) 
and lower (by <30%) N RGB / Nue-ci and N$gb I N B e-ci ratios. The isochrone with A os =0.1 
from PEL dataset reasonably reproduces all the population ratios. Only the Nms/N^sgb+rgb) 
ratio turns out to be ~15% lower than the observed one. 

We conclude that the best agreement with observations (both in terms of evolutionary se- 
quence morphology and star counts) has been obtained by using PEL models computed with 
a mild overshooting (A os =0.1) and r=1.9 Gyr. Also, the required values of distance modulus 
and reddening are fully consistent with those generally adopted for NGC 1978. In order to 



5 Note that the bluest portion of the SGB can be affected by blending. To check this effect, we also defined 
a second box sampling the SGB population (SGB S ), by excluding the bluest region at (V555 — Isu) <0.7. 
The population ratios obtained by using this selection box (and reported in Table 2) are fully consistent with 
the results by using the standard SGB box, suggesting that blending effects (if any) in the SGB population 
have a negligible impact on the results. 
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estimate the overall age uncertainty, we took into account the major error source, namely 
the distance modulus. Hence, we have repeated the best-fitting procedure by using the PEL 
isochrones with mild-overshooting, and varying the distance modulus by ±0.05 and ±0.1 
mag with respect to the reference value of 18.5. A variation of ±0.05 mag still allows a good 
fit of the CMD features with isochrones within =p0.1 Gyr from the reference value of 1.9 Gyr. 
A variation of ±0.1 mag in the distance modulus, does not allows to simultaneously fit the 
He-Clump magnitude level and the extension of SGB, whatever age is selected. Hence, we 
can assign a formal error of ±0.1 Gyr to our age estimate. 



5. Discussion and Conclusions 



The photometric analysis of the ACS-WFC CMD of NGC 1978 presented here pro- 
vided three major results, that can be summarized as follows: i) the firm detection of the 
RGB Bump at V5 55 =19.10±0.10, ii) a new, independent estimate of the cluster ellipticity 
(e=0.30±0.02) and Hi) an accurate measure of the cluster age (r=1.9±0.1 Gyr). 

The detection presented here is the first clearcut detection of the RGB Bump in an 
intermediate age cluster and it confirms the theoretical expectation that this feature also 
occurs in relatively massive stars (the estimated TO mass for this cluster is ~1.5 M , see 
Tab. 1). Note that this result opens the possibility to study the behaviour of the RGB Bump 
as a function of the cluster age and it provides crucial insight on the internal structure of 
intermediate mass stars. 

The high-ellipticity (e=0.30±0.02) of NGC 1978 poses two major questions: i) why the 
LMC clusters in general, and NGC 1978 in particular, are, in average, more elliptical than 
those in the Mil ky Way? ii) why NGC 1978 is more elliptical than the other LMC clusters? 
Goodwin! (119971 ) suggests that the relatively small LMC tidal field can preserve the pristine 
triaxial structure of the clusters, while the strong tidal field of our Galaxy tend to destroy 
it, thus removing at least part of the ellipticity. 

In order to explain the especially high ellipticity of NGC 1978 three main hypothesis have 
been proposed in the past: a merging episode, a rotation effect and an anisotropic velocity 
dispersion tensor (IFischer et al.lll992l ). The mergin g scenario has been proposed because the 
broad RGB from ground-based BVRI photometry (lAlcaino et al.lll999l ) and because the pre- 
limin ary evidence of a metallicity dispersion from high-resolution spectroscopy of two RGB 
stars (IHill et al.ll2000l ) (<5[Fe/H]~0.2-0.3 dex). However, the tiny RGB sequence presented in 
this work as well as the r ecent iron abundance estimate from high-resolution spectra of eleven 
RGB stars presented by iFerraro et al.l (120061 ) definitely excluded any significant metallicity 
spread within the cluster. 
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Our new age estimate (r=1.9±0.1 Gyr) of NG C 1978, coupl e d with the new iron abun- 
dance determination ([Fe/H]=-0.38±0.02 dex) by iFerraro et al.l (120061 ). provide new coor- 
dinates for this cluster in the age-metallicity plane. This is especially important, since the 
correct shape of the AMR in the LMC is still matter of debate: in particular the origin of 
the observed bimodality in the LMC cluster age d istribution has been interpr eted as the 
evidence for two major episodes of star formation. iPagel fc Tautvaisiend (119981 ) computed 
two different AMR semi-empirical models for the LMC, with a continuous star formation 
and with two burst episodes occurred ~14 and 3 Gyr ago, respectively. Fig. 8 shows the 
results of these theoretical predictions. For comparison, the position of N GC 1978 in the 
age- metall ic ity diagram based on i) old metallici ty (— 0.6 < [Fe/H] < — 0.4 lOlszewski et al. 



see 



1984; 



991) ; Ide Freitas Pacheco. Barbuy fc Idiart 



Geisler et al. 



1997; 



Girardi et al 



see 



1998|)) and age (2-3.3 Gyr [Olszewski 



19951) estim ates (grey box), it) the iron abundance 
by Hill et al. ( 2000 ) and age by Bomans et al. ( 19951 ) (open square), and finally in) the most 
recent metallicity by Ferraro et al. ( 20061 ) and age from the present work (big black dot) are 
shown. It is interesting to note that old coordina tes (grey box) ba rely fit with the bursting 
scenario, while the more recent measurements by iHill et al.l (120001 ) place the cluster far be- 
low any model. Our new coordinates are somewhat consistent with both the proposed star 
formation scenarios. 

Similar accurate ([Fe/H], r) coordinates for a significant number of LMC clusters with dif- 
ferent ages and metallicities are urgently needed to disentagle different formation scenarios. 
This is the aim of our ongoing global project. By combining detailed chemical abundance 
(from high-resolution spectra) and ages (from high quality photometry) to a number of pillar 
LMC clusters, we plan to calibrate a suitable age and metallicity scale for the entire LMC 
globular cluster system, with the ultimate goal of providing a robust AMR. 



We warmly thank the anonymous referee for his/her suggestions. This research was 
supported by the Agenzia Spaziale Italiana (ASI) and the Ministero dell'Istruzione, del- 
l'Universita e della Ricerca. 
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Table 1. Age, distance modulus, reddening , Turn-Off mass and magnitude level of the 
RGB-Bump from best-fit BaSTI, PEL and Padua isochrones. 





BaSTI 


BaSTI 


PEL 


PEL 


PEL 


PADUA 




A os =0 


A os =0.2 


A os =0 


A os =0.1 


A OS =0.25 


A OS =0.25 


Age (Gyr) 


1.9 


3.2 


1.7 


1.9 


2.5 


2.2 


(m - M) 


18.47 


18.43 


18.50 


18.50 


18.50 


18.38 


E(B-V) 


0.09 


0.09 


0.09 


0.09 


0.09 


0.07 


M TO (M ) 


1.47 


1.45 


1.49 


1.49 


1.44 


1.45 




19.10 


19.22 


18.73 


18.88 


19.39 


19.44 



Table 2. Theoretical population ratios from BaSTI, PEL and Padua best-fit isochrones, 
and corresponding observed ratios for NGC 1978. 



Population Ratio 


BaSTI 


BaSTI 


PEL 


PEL 


PEL 


PADUA 


Observed 




A os =0 


A os =0.2 


A os =0 


A os =0.1 


A OS =0.25 


A OS =0.25 




MS/(RGB+SGB) 


2.40 


6.05 


2.10 


3.27 


7.06 


6.92 


4.00±0.40 


SGB/He-Cl 


4.16 


1.63 


4.90 


2.01 


1.62 


1.05 


2.03±0.22 


RGB/He-Cl 


1.95 


1.07 


2.41 


1.64 


1.03 


0.78 


1.45±0.19 


RGB/SGB 


0.58 


0.65 


0.49 


0.81 


0.64 


0.75 


0.71±0.10 


SGB a /He-Cl 


1.39 


0.64 


1.43 


0.85 


0.61 


0.36 


1.12±0.12 


MS/(RGB+SGB a ) 


4.14 


9.61 


4.00 


4.75 


11.43 


11.29 


5.43±0.42 
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Fig. 1.— (V, V-I) CMD of the LMC cluster NGC 1978, obtained with ACS@HST (only first 
chip). The dashed line indicates the saturation level. 
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Fig. 2.— (V, V-I) CMD of the outer region (r>140" from the cluster center) of NGC 1978, 
as obtained with ACS@HST (only second chip). 
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Fig. 3. — Completeness curves computed in two radial sub-regions of NGC 1978. The black 
points indicate the value of the = parameter calculated for each magnitude bin. The 

IV sim 

completeness curves have been computed for V555 <17.6, corresponding to the saturation 
level. Tipycal errorbars are also indicated. 
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Fig. 4. — Left: the bright portion of the CMD of NGC 1978 (grey points) with the selected 
RGB stars (black points). The arrow indicates the magnitude of the RGB Bump. Right: 
differential (upper panel) and integrated (lower panel) LFs, computed for the RGB stars, 
excluding the He-Clump and AGB populations. The arrow in the upper panel indicates the 
position of the RGB Bump. The dashed lines in the lower panel are the linear fit to the 
regions above and below the RGB Bump. 
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Fig. 5. — Upper panel: the map of NGC 1978 with the isodensity contours; central panel: 
the best fit ellipses of the isodensity contours; lower panel: ellipticity of the best fit ellipses 
as a function of the semi-major axis in arcsecond. The horizontal dashed line indicate the 
mean value. 
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Fig. 6. — Best-fit theoretical isochrones overplotted to the observed CMD of NGC 1978 (only 
stars at r>20" from the cluster center are plotted) obtained with theoretical isochrones: each 
panel shows a different model and the corresponding A os value and age. Reddening and 
distance moduli for each model are reported in Table I. 
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Fig. 7.— The bright portion of the NGC 1978 CMD (only stars at r>20" from the cluster 
center are plotted) with the selection boxes adopted to sample the MS, SGB, RGB and 
He-Clump populations. 
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Fig. 8. — Theoretical predictions for the LMC AMR computed by iPagel fc Tautvaisiene 
(119981 ): continuos line refers to a AMR obtained assuming a continuous SF, and dashed line 
corresponding to a bursti ng model. The blac k point indicates the position of NGC 1978 
using the metallicity from lFerraro et al.l ((20061) and the age de rived in th i s stud y. The open 
square refers the po sition for this cluster using metallicity by iHill et al.l (120001 ) and age by 



Bomans et al.l (119951 ). The grey box shows previous metallicity (—0.6 < [Fe/H] < —0.4) and 



age (2-3.3 Gyr) estimates from low resolution data. 



